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ABSTRACT 



We present the attempt to study the problem of the estimates of higher-order perturbative 
corrections to physical quantities in the Euclidean region. Our considerations are based on the 
application of the scheme-invariant methods, namely the principle of minimal sensitivity and 
the effective charges approach. We emphasize, that in order to obtain the concrete results for 
the physical quantities in the Minkowskian region the results of application of this formalism 
should be supplemented by the explicite calculations of the effects of the analytical continuation. 
We present the estimates of the order 0(a^) QCD corrections to the Euclidean quantities: the 
e + e~-annihilation D-function and the deep inelastic scattering sum rules, namely the non- 
polarized and polarized Bjorken sum rules and to the Gross-Llewellyn Smith sum rule. The 
results for the D-function are further applied to estimate the 0(aj) QCD corrections to the 
Minkowskian quantities R(s) = <7 io t(e + e~ — > hadrons)/(x(e + e~ — > fi + ^~) and R T = T(r — > 
v T + hadrons)/r(r — > u T V e e). The problem of the fixation of the uncertainties due to the O(af) 
corrections to the considered quantities is also discussed. 
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I. During the last few years, essential progress has been achieved in the area of the calcu- 
lation of the next-next-to-leading order (NNLO) QCD corrections to the number of physical 
quantities. Indeed, the complete NNLO O(a^) QCD corrections are known at present for the 
characteristics of e + e~ — > hadrons process 0, r — ► z/ r +hadrons decay [[| Z° — ► hadrons pro- 
cess Q and for the deep inelastic scattering sum rules, namely the non-polarized Bjorken sum 
rule (BjnSR) ||, the Gross-Llewellyn Smith sum rule (GLSSR) and the polarized Bjorken sum 
rule (BjpSR) |J. Amongst the physical information provided by these results is the estimate 
of the theoretical uncertainties of the corresponding next-to-leading-order (NLO) perturbative 
QCD predictions for R(s) = a tot (e + e~ — > hadrons)/cx(e + e~ — > 0, T(Z° — > hadrons) 

(see. e.g., §), R T = T(r -> v T + hadrons) /T(r -> z/ T z7 e e) §, BjnSR 0, GLSSR and BjpSR 
B 

After gaining from the understanding of the phenomenological meaning of the effects of the 
NNLO corrections to observable quantities one can of course stop at this level. However, in 
view of the fact that the precision of the experimental data for R(s), R T , T(Z° — > hadrons) 
and the DIS sum rules is continuously increasing, the questions about the possible values of 
the effects of non-calculated higher order terms is frequently raised both by theoreticians and 
experimentalists. This information can be useful in the studies of several important problems. It 
allows to present the concrete numerical estimate of the theoretical uncertainties of the NNLO 
QCD approximations. Moreover, any new information, even empirical, about the possible 
values of the higher-order terms might be useful in understanding of the general structure of 
the perturbative expansions for different quantities in the concrete renormalization schemes 
and different energy regions. 

A certain step in the direction of the estimates of the order 0(af) QCD corrections was 
made in the case of R T in Ref. 0. These estimates are based on the tendency, observed in Ref. 



ll|], of the scheme-dependent uncertainties of the perturbative QCD predictions for R(s) and 
R T to decrease as a result of taking into account the order 0(a;?)-ternis. This effect already 
occurs at the 0(a 2 s ) level. The inclusion of the 0{a z s ) corrections fill confirms this effect and 



makes it more vivid. The foundations of Ref. |TT[ were further improved in the process of 
phenomenological studies of the QCD predictions for R T , using the explicite separations of the 
expressions for the n-th order coefficients r T n to the n-th order coefficients d n of the Euclidean 
.D-function and the certain contributions, which can be explicitely taken into account after 
completing the calculations in the k < n — 1 order of the perturbation theory [ ]E| . 

Notice, that starting from the NNLO level these contributions contain 7r 2 -factors which 
appear as the result of analytic continuation to the physical region. As was demonstarted in 



Refs. |L3|, |l^|,p2| in order to achieve better convergence of the corresponding approximants it 
is essential to treat them seperately from the coefficients of the D-function. Even the seperate 
renormalization-group ]15] inspired summation of all additional terms in the expression for R(s) 



[|TJ] and R T [|T3|, [T2j was discussed. 

This work is the attempt to address the important questions of the study of the higher- 
order perturbative approximations to physical quantities. We will use the approach based on 
the imrovement formula |L6| of re-expansion of the expressions for the quantities obtained in the 
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N-th order of perturbation theory within the principle of minimal sensitivity (PMS) [TfJ and the 
effective charges (ECH) approach |T7j , which is equivalent a posteriori to the scheme-invariant 



perturbation theory |fL8| . The main aim of this work, which is based on the considerations of 



Refs. [19, 20], is to present the observations obtained in the process of the study of the problem 
of the possibilities of estimating higher-order QCD corrections to R(s), R T and deep-inelastic 
scattering sum rules using the scheme-inariant procedures. 

Of course, this approach to the estimation of the values of the uncalculated higher- order 
terms cannot be considered as the alternative to the direct analytical or numerical calculations. 
However, we hope that this method can give the concrete prescription for impression about 
the possible theoretical errors due to variation of the higher order terms. The first message 
came from the fact that the rather bold-guess application of this method for the estimate of the 
four-loop corrections to the expression for (g — 2) e pl| gave results, which turned out to be in 



surprisingly good agreement with the latest results of the direct numerical calculations of Ref. 
2lJ 0. The second argument in favour of this procedure came from its successful application 



for the analysis of the Drell-Yan cross section at the 0(al)-leve\ [2~4j]. Moreover, as will be 



demonstrated in our work, the re-expansion formalism of Ref. ]T6| is also working quite well in 
QCD at the NNLO for at least three independent quantities, namely for the e + e~-annihilation 
.D-function and the non-polarized and polarized Bjorken sum rules. Therefore, we prefer to 
consider all these facts as an arguments in favour of more detailed studies of the intrinsic 
features of this approach. 



In this work we further develop the re-expansion formalism of Ref. [IB], deriving new terms 
in the corresponding improvement formula. The previously-known terms are used to obtain 
the estimates of the next-to-next-to-next-to-leading order (N 3 LO) QCD corrections to the D- 
function, the non-polarized Bjorken sum rule and polarized Bjorken sum rule, which is closely 
related by the structure of the corresponding perturbative series to the Gross-Llewellyn Smith 
sum rule. We will use the results obtained for the D-function to estimate the effects of the 
N 3 LO corrections to the perturbative series for the Minkowskian quantities R(s) and R T by 



adding the explicitly calculable terms, previously discussed in Refs. [25, 12|, supplementing thus 



the related considerations of Refs. [16], ]£6) by the additional input information. The derived 
new term in the improvement formula is applied to touch on the problem of fixing the values of 
the 0(af)-corrections to the analysed quantities. We also specify the scheme which is related 
to our considerations and emphasize that the structure of the perturbative expansions in the 
Euclidean and Minkowskian regions differ essentially. 

II. Consider first the order 0(a N ) approximation of a Euclidean physical quantity 

N-l 

D N = d a(l + dia 1 ) (1) 

i=l 

with a = a s /n being the solution of the corresponding renormalization group equation for the 

4 A detailed re-consideration of the (g — 2) M analysis of Ref. Q| and its generalization to the five-loop order 
will be presented elsewhere p3[. 
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/3-function which is defined as 



da 
<9/i 2 



N-l 

(3(a) = -f3 Q a 2 (l + ^ c i a% 

i=l 



(2) 



In QCD, in the process of the concrete calculations of the coefficients di, i > 1 and q, % > 2, the 
MS* scheme is commonly used. However, this scheme is not the unique prescription for fixing 
the RS ambiguities (for the recent discussions see e.g., Ref. |26|] ). 



The PMS |16| and ECH jl7| prescriptions stand out from various methods of treating 
scheme- dependence ambiguities. Indeed, they are based on the conceptions of the scheme- 
invariant quantities, which are defined as the combinations of the scheme- dependent coefficients 
in Eqs. (1) and (2). Both these methods pretend to be "optimal" prescriptions, in the sense 
that they might provide better convergence of the corresponding approximations in the non- 
asymptotic regime, and thus allow an estimation of the uncertainties of the perturbative series 
in the definite order of perturbation theory. Therefore, applying these "optimal" methods, one 
can try to estimate the effects of the order 0(a Ar+1 )-corrections starting from the approxima- 
tions D < ^ t (a opt ) calculated in a certain "optimal" approach 



2T| , P7| . This idea is closely 



related to the QED technique of Ref. [p8| |, which was used to predict the renormalization-group 
controllable Zn(m^/m e )-terms in the series for (g — 2) M from the expression of (g — 2) M through 
the effective coupling constant a(m^/m e ). In our work work we are using this technique to 
estimate also the constant terms of the higher-order corrections in QCD. 



Let us following the considerations of Ref. |16| re-expand D c ^ t (a op t) in terms of the coupling 



constant a of the particular scheme 



D 



opt I 
N \ 



l opt) 



D N (a)+5D 



°pt „ N+1 



N 



a 



(3) 



where 

5Df = n^d,, a) - n N (df, cf) (4) 

are the numbers which simulate the coefficients of the order 0(a Ar+1 )-corrections to the physical 
quantity, calculated in the particular initial scheme, say the MS'-scheme. The coefficients 
can be obtained from the following system of equations: 





dr 



N 



n N a N+1 



d 



^-(D N + Q N a N+1 ) 

OCi 



0(a N+2 ), i>2 



(5) 



where the parameter r = f3 Q £n([i 2 /A 2 ) represents the freedom in the choice of the renormal- 
ization point fi. The conventional scale parameter A will not explicitly appear in all our final 
formulas. 

The explicit form of the coefficients f2 2 and ^ in which we will be interested can be obtained 
by the solution of the system of equations (5), following the lines of Ref. 
the final already known expressions 



16]. We present here 
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il 2 = rf rfl(Cl + G?i), 



(6) 
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1 

£7 3 = d di(c 2 - -cidi - 2d\ + 3d 2 ) . (7) 
and the new term f2 4 evaluated by us: 

fi 4 = ^(3c 3 rfi + c 2 d 2 - Ac 2 d\ + 2c x d x d 2 - Cl d 3 + \Ad\ - 28d 2 d 2 + 5d 2 2 + 12did 3 ) (8) 
which reproduces the renormalization-group controllable logarithmic terms at the five-loop level 



It should be stressed that in the ECH approach df CH = for alH > 2. Therefore one gets 
the following expressions for the NNLO and N 3 LO corrections in Eq. (3): 

6D% CH = n 2 (d 1 ,c L ) (9) 

5Df CH = Q 3 (d 1 ,d 2 ,c 1 ,c 2 ) (10) 

SDf CH = fi 4 (di, da, d 3 , c 1; c 2 , c 3 ) (11) 
where Q 2 , f2 3 and Q 4 are defined in Eqs. (6), (7) and (8) respectively. 

It is worth emphasizing that the general expressions for the correction coefficients 5D^ = 
d d N in Eqs. (9), (10), (11) can be obtained from the following exact relation for the process 
dependent but scheme-independent quantities 

° N - ° N +d N -^ (12) 



N-l N-l JV d, 



o 



which are related to the coefficients c% CH of the ECH /^-functions defined as P e ff(o,ECH) = 
— Po^ECH(^-~^ c i a ECH + J2i>2 c f CHa ECH)- Therefore, the expressions for the corrections SD§ CH 
are the exact numbers which are related to the special scheme. This scheme is identical to the 
MS scheme at the lower order levels and is defined by the condition cn = cff CH at the N-th 
order, where cfj CH are considered as unknown numbers. 

In order to find similar corrections to Eq. (3) in the N-th order of perturbation theory starting 
from the PMS approach it is necessary to use the relations obtained in Ref. |30| between 



the coefficients df MS and cf MS (i > 1) in the expression for the order 0(ap MS ) approximation 
Dff MS (apMs) of the physical quantity under consideration. The corresponding corrections have 
the following form: 

5D™ 13 = 5D* CH + ^ (13) 

8Di MS = 5Di CH . (14) 
Notice the identical coincidence of the N 3 LO corrections obtained starting from both the PMS 



and ECH approaches. A similar observation was made in Ref. |2T] using different (but related) 
considerations. 

In the fourth order of perturbation theory the additional contribution to 5D^ MS has more 
complicated structure. Indeed, the expression for VL i (df MS ) cf MS ) in Eq. (4) reads: 

O (d PMS r PMS \ - r PMS 4 ( r PMS\2 5 2 PMS , 7 4n /-, r\ 

"4l«i ,c { ) - yi^ c i c 3 ~ 81^ 2 > ~ 81 1 2 648 
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where 

c P 2 MS = \(d 2 + c 2 -d\- ah + ^c\) + 0(a) (16) 

and 

cf M5 = 4(d 3 + lc 3 - c 2 d x - 3d ± d 2 + 2d\) + -ci(d 2 + c 2 + 3d 2 - c x d x + — c?) + 0(a). (17) 
2 2 108 

The expressions for Eqs. (15)- (17) are pure numbers, which do not depend on the choice of 

the initial scheme. Note that we have checked that in the case of the consideration of the 

perturbative series for (g — 2) M the numerical values of Q^(df MS ,cf MS ) are small and thus 

the a posteriori approximate equivalence of the ECH and PMS approaches is preserved for the 



quantities under consideration at this level also p3fl . We think that this feature is also true in 
QCD. 

III. Consider now the familiar characteristic of the e + e _ — > 7 — ► hadrons process, namely 
the Z)-function defined in the Euclidean region: 

R(s) 



Its perturbative expansion has the following form: 

D(Q 2 ) = 3SQj[l + a + J2 dia i+1 ] + (ZQ f ) 2 [d 2 a 3 + 0(a 4 )] (19) 



i>l 

\2 



where Qf are the quark charges, and the structure proportional to (SQ/) comes from the 
light-by-light-type diagrams. The coefficients di and d 2 , d 2 were calculated in the MS'-scheme 
in Refs. @ and |], respectively. They have the following numerical form: 

dfS w 1.986 - 0.115/ 

dfs Rj 18.244 - 4.216/ + 0.086/ 2 , d 2 w -1.240 . (20) 



Following the proposals of Ref. [5I|, we will treat the light-by-light-type term in Eq. (19) 



separately from the "main" structure of the D-function, which is proportional to the quark- 
parton expression D QP (Q 2 ) = 3T.Q 2 . In fact, one can hardly expect that it is possible to 
"predict" higher-order coefficients d^i > 3 of the second structure in Eq. (19) using the 
only explicitly-known term d 2 . Therefore we will neglect the light-by-light-type structure as 
a whole in all our further considerations. This approximation is supported by the relatively 
tiny contribution of the second structure of Eq. (19) to the final NNLO correction to the 
D-function. 

The next important ingredient of our analysis is the QCD /5-function (2), which is known 
in the MS-like schemes at the NNLO level |32fl . Its corresponding coefficients read 



Po = (11 ~\f)\~ 2 -75 -0.167/ 



153-19/ 

[D99 f + 395 f 2 

(21) 



66-4/ 

m 77139 - 15099/ + 325/ 2 



9504 - 576/ 
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Using now the perturbative expression for the -D-function, one can obtain the perturbtive 
expression for R(s), namely 



R(s) = 3SgJ[l + a s + na l s +1 ] + (EQ f ) 2 [f 2 a 3 s + 



(22) 



where a s = a s /n, and 



ri = di 

r - d ^ f - d 
r 2 — a 2 — , r 2 — a 2 



r 3 = d 3 -Tc f3 (d 1 + -ci) . 

o 



(23) 



The corresponding 7r 2 -terms come from the analytic continuation of the Euclidean result for 
the D-function to the physical region. The effects of the higher-order 7r 2 -terms were discussed 
in detail in Ref. |25[|. For example, the corresponding espression for the r4-term reads : 



7 1 7T^ 

r 4 = d 4 - 7r 2 /3 2 (2d 2 + -cidi + -c\ + c 2 ) + y^ 4 



(24) 



The perturbative expression for R T is defined as 

rM? ds 



Rr = 2j o w (l- s/M*Y (1 + 2s/M')R(s) ~ 3[1 + a T + £ r[< +1 ] 



(25) 



where a T = a s (M 2 )/7r and i?(s) is i?(s) with with / = 3, (SQ/) 2 = 0,3S<5 2 substituted for 



3S | | 2 and | V ud | 2 + | V t 



1. 



It was shown in Ref. 1T3] that it is convenient to express the coefficients of the series (25) 
through those ones of the series (19) for the D-function in the following form : 



d? s {f = 3)+gi{f = 3) 



d% IS (f = 3)+g 2 (f = 3) 



di vrs (f = 3)+9 3 (f = 3) 



(26) 



where in our notations 



9i = -Poh = {f = 3) = 3.563 

g 2 = -{2d l + c l ]f3 I l +f3ll 2 = (f = 3) = 19.99 

93 = -[3d 2 + 2d 1 c 1 + c 2 }p I 1 + [3d 1 + ^c 1 }P^I 2 -Pll 3 = (f = 3) = 78.00 



(27) 



and Ik are defined and calculated in Ref. [O. Their analytical expressions read: 1\ = —19/12 
I 2 = 265/72 - tt 2 /3 and J 3 = -3355/288 + 19vr 2 /12. One of the pleasant features of Eqs 
is that they are absorbing all effects of the analytical continuation. 
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Following the lines of Ref.|T2|| we derive the corresponding expression for the coefficient r\: 

rl = df^(f = 3) + g,(f = 3) (28) 

where 

3 

g A = -[Ad 3 + 3d 2 ci + 2d lC2 + c 3 ]p Ii + [6d 2 + lc x dx + -c\ + 3c 2 ]^2 
13 

-[4 C / 1 + y c 1 H 3 J 3 + /5 4 J 4 = (/ = 3) 

= 3.562c 3 (/ = 3) + 14.247rf 3 (/ = 3) - 466.73 (29) 
and I A = 41041/864 - 265tt 2 /36 + tt 4 /5 « -5.668. 

In order to estimate the values of the order 0(a 3 ), 0(a 4 ) and 0(a 5 ) corrections to R(s) 
and R T , we will apply Eqs. (6) - (11) in the Euclidean region to the perturbative series for the 
/^-function and then obtain the estimates we are interested in using Eqs. (23), (24), (26)-(29). 
This is the new ingredient of this analysis, which distinguishes it from the related ones of Refs. 

We now recall the perturbative expression for the non-polarized Bjorken deep-inelastic scat- 
tering sum rule 



BjnSR = C F! p ' vp (x, Q 2 )dx = 1 - -ail + V d l( j) (30) 
Jo 3 f£ 



where the coefficients d\ and d% are known in the MS scheme from the results of calculations 
of Ref. and Ref. || respectively: 

dfs ~ 5.75 _ 0.444/ (31) 

dfs « 54.232 - 9.497/ + 0.239/ 2 . (32) 
The expression for the polarized Bjorken sum rule BjpSR has the following form: 

BjpSR = C gT~ en (x, Q 2 )dx = - | — | [1 - a(l + V d^)} (33) 

3 g v {>1 



where the coefficients di and d 2 were explicitly calculated in the MS scheme in Refs. [H| and 
respectively. The results of these calculations read 



df s w 4.583 - 0.333/ (34) 



df b w 41.440 - 7.607/ + 0.177/' • (35) 

It should be stressed that since deep inelastic scattering sum rules are defined in the Euclidean 
region, we can directly apply to them the methods discussed in Section 2 without any additional 
modifications. 
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It is also worth emphasizing that, in spite of the identical coincidence of the NLO correction 
to the Gross-Llewellyn Smith sum rule GLSSR = (1/2) f F^ p+vp (x, Q 2 )dx with the result of 
Eq. (34) jlOfl , the corresponding NNLO correction differs from the result of Eq. (35) by 
the contributions of the light-by-light-type terms typical of the GLSSR ||. Since these light- 
by-light-type terms appear for the first time at the NNLO, it is impossible to predict the 
value of the light-by-light-type contribution at the N 3 LO level using the corresponding NNLO 
terms as the input information. However, noticing that at the NNLO level the corresponding 
light-by-light-type contributions are small M, we will assume that the similar contributions are 
small at the N 3 LO level also. Only after this assumption can our estimates of the NNLO 
and N 3 LO corrections to the BjpSR be considered also as the estimates of the corresponding 
corrections in the perturbative series for the GLSSR . Note, that the Pade predictions of the 
order 0(a 4 ) contributions to the GLSSR [Q, which do not take into account the necessity of 
the careful considerations of the light-by-light-type terms, deserve more detailed considerations. 
To our point of view, it is better to neglect the small light-by-light contribution as the whole 
and to consider the problem of estimates of the higher order corrections to BjpSR and GLSSR 
simultaniously, but not seperately, as was done in Ref. p3| . 



IV. The estimates of the coefficients of the order 0(a 3 ) and 0(a A ) QCD corrections 
to the D-function, R(s), BjnSR and BjpSR/GLSSR, obtained following the discussions of 
Section 2 with the help of the results summarized in Section 3, are presented in Tables 1-4 
respectively. Due to the complicated /-dependence of the coefficients ^2,^3 in Eqs. (6) and 
(7), we are unable to predict the explicit /-dependence of the corresponding coefficients in the 
form respected by perturbation theory. The results are presented for the fixed number of quark 
flavours 1 < / < 6. The estimates of the NNLO corrections, obtained starting from both 
the ECH and PMS approaches, are in qualitative agreement with the results of the explicit 
calculations. The best agreement is achieved for / = 3 numbers of flavours. 

Using the results of Table 1 and Eqs. (25), (26) we get the following estimates of the NNLO 
coefficients R T : (r^^Qfj ~ 25.6 and {r T 2 )p MS « 24.8 . One can see the agreement with the 
explicitly calculated result (r^jjg = 26.366 . Considering this agreement as the additional a 
posteriori support of the methods used, we use the estimate of the N 3 LO coefficient for the 
/^-function with / = 3 numbers of flavours as presented in Table 1, namely: 

df (/ = 3) = 27.5 (36) 

and estimate the value of the N 3 LO coefficient of R T : {r^) est « 105.5 . 

As follows from Eq. (12) the results given in Tables 1-4 correspond to the ECH-inspired 
variant of the MS scheme. However, we hope that there is some meaning in the relation = 
c^ CH used to obtain these numbers. Indeed, it is known from the explicite QED calculations of 
Refs. []34] that the difference between the numerical values of the higher order coefficients of the 



QED /5-function in different schemes is decreasing at the four-loop level. So, we do not exclude 
the situation that the difference between scheme-invariants c^ CH of different QCD quantities, 
which probably have the factorial growth cfj CH ~ AN\, will also decrease in the higher orders 
of perturbation theory, and that it might happen that cn ~ c% CH . If so, the estimates of Tables 
1-4 reveal the structure of the perturbative series for physical quantities in the fixed scheme. 
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In order to address the problem of fixing the values of the 0(a 5 ) QCD corrections to the 
considered Euclidean quantities we apply Eqs. (8), (11) with the explicitely calculated NLO 
and NNLO coefficients di, d 2 , ci, c 2 in the MS scheme and use the determined from Eqs. (7), 
(10) estimates of the N 3 LO coefficients d 3 . In order to obtain the estimates of the next-to- 
next-to-next-to-next- to-leading order (N 4 LO) coefficients of R(s) and R T , which are related to 
the N 4 LO coefficient d\ of the D-function, the explicitly calculated terms in Eq.(24) and Eq. 
(29) are taken into account. 

However, the expressions for f2 4 in Eq. (8) and g 4 in Eq. (29) depend also on the four-loop 
coefficient c 3 of the QCD /^-function, which is unknown at present. Therefore, in Tables 1-4 
we present the estimates of the combinations g? 4 — <iic 3 and r 4 — r\C^. The existing uncertainty 
in the value of N 3 LO coefficients <i 3 is fixed by the assumption that the real values of these 
coefficients do not significantly differ from the N 3 LO estimates, obtained by us. However, in 
the case of the D-function with / = 3 numbers of flavours we present also the more detailed 
expression of <i 4 , which follows from Eqs. (11), (8): 

(d±)Tc H (f = 3) = l-64c 3 (/ = 3) + 5.97d 3 (/ = 3) - 52.8. (37) 

Taking into account the negative contributions into the coefficient g 4 in Eq. (28) we get the 
following estimate of the N 4 LO coefficient of R T : 



1. 



f CH = 5.2c 3 (/ = 3) + 20.22rf 3 (/ = 3) - 519.5 . (38) 



Of course, the best way of fixing the value of c 3 would be its explicite calculation. However, at 
the current level of art one is free to invent his own way of fixing this part of the ambiguities 
of the less substantiated than N 3 LO estimates of the N 4 LO terms []. We will use here the 
bold guess estimate c 3 (/ = 3) ~ c 2 {f = 3) 2 /ci(/ = 3) ~ 11 , which is motivated by the good 
agreement of our results of Eq. (36) with the approximation d 3 « d\jd\ « 25, previously used 
in Ref. (12| to fix the value of the N 3 LO coefficient of the Z)-function for / = 3 . Combining 



it with the estimate of Eq. (36) we get the following estimate of the N LO coefficients of R 



{r T A ) est w 93. 



V. We are now ready to discuss the main otcomes of our analysis. 

1. The estimates of the NNLO corrections to the D-function, BjnSR and BjpSR/GLSSR 
are in qualitative agreement with the results of the explicit calculations of Refs. |], ]2|, || 
and and respect the tendency of the corresponding coefficients to decrease with increasing 
number of flavours. 

2. The best agreement of the NNLO estimates with the exact results is obtained for the 
case / = 3. This fact supports the application of the method used for estimating the NNLO 
and N 3 LO corrections to R T . 

3. Notice that since the methods used correctly reproduce the renormalization-group- 



controllable terms [2S|, [IB], the transformation of the coefficients di, d 2 from the MS scheme 

5 Another kind of the uncertainties is related to the possible deviations of the real values of the N 3 LO 
corrections from the presented by us estimates. 
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to other variants of the MS-like scheme will not spoil the qualitative agreement with the results 
of the explicit calculations. 

4. The comparison of the results of Table 1 with those of Table 2 demonstrate that the 
7r 2 effects give dominating contributions to the coefficients of R(s). 

5. Amongst other outcomes are the estimates of the N 3 LO and N 4 LO corrections to R T 
and R(s) for / = 5 numbers of flavours. Taking a s [Mz) ~ 0.12, we get the estimate of the 
corresponding N 3 LO contribution to both T(Z° — ► hadrons) and T(Z° — > qq): 

($T z o) N 3 LO « -97(a(M z )) 4 « -2 x 10~ 4 . (39) 

It is of the order of magnitude of other corrections included in the current analysis of LEP 
data (see, e.g., fl35| ). Using the estimates presented in Table 2 and assuming that c 3 (/ = 5) = 
c 2 (/ = 5) 2 /ci(/ = 5) ~ 1.715 one can also estimate the corresponding N 4 LO contributions, 
namely 

(ST z o) N * LO « 70(a(M z )) 5 « 5 x lO" 6 (40) 

The way of fixation of the value of c 3 used above is not applicable for the case of / = 6 since 
we expect that in this case its real value is negative. 

6. Taking a s (M T ) pa 0.36 ||, we get the numerical estimate of the N 3 LO contribution to 

R T \ 

(5R T ) N s LO ^ 105.5 a\ pa 1.8 x 10~ 2 . (41) 
It is larger than the recently calculated power-suppressed perturbative [3B| and non-perturbative 



3^] contributions to R T . The possible contributions of the N A LO corrections are significantly 
smaller: 

(6R T ) N 4 LO ^ 93a 5 T pa 1.8 x 10~ 3 . (42) 

7. Notice, that in spite of the possible N\ growth of the coefficients of the D-function, 
the additianal 7r 2 -dependent contributions are shadowing down this possible growth in the 
Minkowskian region. In order to consider this effect in higher order levels we applied the RG- 
technique and derived the corresponding relations between the order 0(a 6 ) coefficients of the 
/^-function and R(s) and R T . In the case of R(s) the relation reads 

2 27 y 7 

r 5 = d 5 - y/3o(10c/ 3 + —c x d 2 + 4,c\di + -cic 2 + 8c 2 di + -c 3 ) 

+ ^ 4 (5d 1 + ^ Cl ). (43) 

In the case of R T the correction term to the relation = ds(f = 3) + g$(f = 3) is more 
complicated 



g 5 = -[5d 4 + Acid 3 + 3c 2 d 2 + 2c 3 di + C4]f3 h 

27 . ,7 7 



+ [10d 3 + — Cl d 2 + 4c 2 d! + 8^ + -c x c 2 + -c 3 ]/3 2 / 2 

AT Q t\ TT 

-[I0d 2 + —ah + —c\ + Qc 2 }P 3 I 3 + Mi + I ^c 1 ]/3 4 J 4 - (3 5 I 5 (44) 
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where 7 5 = -2479295/10368 + 16775vr 2 /432 - 19tt 4 /12 w -10.113. Keeping in the expression 
for g 5 (f = 3) all explicitely unknown coefficients we get 



g 5 (f = 3) = 17.8d 4 (/ = 3) + 45.07d 3 (/ = 3) + 3.56c 4 (/ = 3) + 18.6c 3 (/ = 3) - 8455. (45) 

Notice the appearence in the expression for g 5 (f = 3) of the huge negative coefficient. The 
derived expressions can mean that the asymptotic structure of the perturbative expansion in 
the Minkowskian region differs essentially from the one in the Euclidean region. 

8. The qualitative agreement of the results of Tables 3 and 4 for BjnSR and BjpSR with 



the corresponding Pade estimates of Ref. [ 33 1 can be considered as the argument in favour of 
the applicability of both theoretical methods in the Euclidean region for the concrete physical 
applications. Let us stress again that in the process of these applications the light-by-light-type 
structures, contributing to the GLSSR, should be treated seperately. 

9. Our estimate for d^f = 3) of Eq. (36) is more definite than the one presented in Ref. 
3^], namely d 3 (/ = 3) = 55^4, and than the bold guess estimate d 3 (/ = 3) = ±25, given 



in Ref. pi. The result is in good agreement with the "geometric progression" assumption of 



Ref. [0. The related estimate of the N 3 LO contribution to R T (see Eq. (41)) is more precise 
than those presented in Refs. H and |12|, namely 5R T = ±130 a 4 [|J and 5R T = (78 ± 25)a 4 



[|T^], and is smaller than the result of applying the Pade resummation technique directly to R T , 
namely 5R T = 133a 4 [0,11. 



10. The application of the Pade resummation technique to R(s) pqj , stimulated by the 
previous similar studies of Ref. [4~0|, gives less definite estimates than the results of applications 
of our methods (compare the estimate 5R(s) = (— 49^4o) ° 4 [§H with the result 5R(s) = —97a 4 
from Table 2 and Eq. (39)). Probably, this fact is related with the problems of applicability 
of the Pade resummation technique to the sign variation perturbative series for R(s) in the 
MS'-scheme. 

11. Let us emphasize again, that the used by us method was applyed directly to the 
Euclidean quantities, namely to the D-function and deep-inelastic scattering sum rules. The 
results for R(s) and R T were obtained from the ones of the D-function after taking into account 
explicitely calculable terms of the analytic continuation to the Minkowskian region. In principle, 
one can try to study the application of the procedure used by us directly to R(s) and R T in 
the Minkowskian region. We checked that in the case of applications of Eqs. (6)-(8) with 
the substitution of dj = or d; t = r\ one can reproduce the same values of the Euclidean 
coefficients dj but the structure of the high order 7r 2 -terms will not agree with the explicitly 
known results. A similar problem will also definitely arise in the case of more rigorous studies 
of the applicability of the Pade resummation methods for the Minkowskian quantities R(s) and 



R T (for the studies existing at present see Refs. |38|, |39|). We think that these analyses should 
be supplemented by the construction of the Pade approximants directly to the D-function. In 
our case the solution of this problem lies in the necessity of modifications of Eqs. (6)-(8) in 
the Minkowskian region by adding concrete calculable 7r 2 -dependent and scheme-independent 
factors. After these modifications it is possible to reproduce our results for R(s) and R T by 
means of application of a similar technique directly in the Minkowskian region. 
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12. Our estimate for d^(f = 3) was recently supported by the phenomenological analysis 
of the ALEPH data for R T jy]]. The above presented considerations of the diffence between 
the structure of the perturbative series for the D-funcion and R(s) already stimulated the 
reconsiderations of the applications of the Pade approximants for the estimates of the higher- 
order coefficients of R(s) and R T |42j| . Note, however, that we are realizing that at the present 
level of understanding both methods suffer from the lack of more rigorous dynamical information 
about the analytical structure of the expanded in the perturbative series functions both in the 
Euclidean and Minkowskian regions. We hope that the considerations presented above will 
stimulate future more detailed studies of the related problems. 

13. The interesting fact is that the presented NNLO estimates are in better agreement with 
the explicite results for / = 3 numbers of flavours than say for / = 5 numbers of flavours. It 
should be stressed again that the basic assumption of our approach is the relation cff H = c^. 
However, it is known that on the contary to the coefficients the scheme-invariants r - Li 11 



-N 



contain additional factors of order / 1 |43j , though with small coefficients. For larger values 



of / these terms become more and more important. This might lead to the additional sourse 
of the violation of the basic relation cfj CH = cn of the re-expansion procedure used by us. 
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f 

J 


df 




( d% s ^} pa/c 


d% st 

6 


(dT 1 ) fcn — c-idi 

V ^4 / Hi O i2 O ^ X 


1 


14.11 


7.54 


7.70 


75 


476 


2 


10.16 


6.57 


7.55 


50 


260 


3 


6.37 


5.61 


6.40 


27.5 


111 


4 


2.76 


4.68 


5.27 


8.4 


22.7 


5 


-0.69 


3.77 


4.16 


-7.7 


-13.2 


6 


-3.96 


2.88 


3.08 


-21 


-2.76 



Table 1 : The results of estimates of the NNLO, N 3 LO and N 4 LO corrections in the series 
for the .D-functions. 



/ 




(r e 2 st )ECH 




r .est 
'3 


(rl st ) E cH ~ c z di 


1 


-7.84 


-14.41 


-14.25 


-166 


-1750 


2 


-9.04 


-12.65 


-11.67 


-147 


-1161 


3 


-10.27 


-11.04 


-10.25 


-128 


-668 


4 


-11.52 


-9.59 


-9 


-112 


-263 


5 


-12.76 


-8.32 


-7.93 


-97 


67 


6 


-14.01 


-7.19 


-6.99 


-83 


330 



Table 2: The results of estimates of the NNLO, N 3 LO and N 4 LO corrections in the series for 
R(s). 



f 


df 


(df^ECH 


(df l )pMS 


d% st 


(dl st )ECH — Csdi 


1 


44.97 


39.62 


40.78 


424 


4127 


2 


36.19 


33.28 


34.26 


303 


2613 


3 


27.89 


27.37 


28.16 


200 


1474 


4 


20.07 


21.91 


22.50 


114 


664 


5 


12.72 


16.91 


17.30 


44 


138 


6 


5.85 


12.39 


12.59 


-10 


-145 



Table 3: The results of estimates of the NNLO, N 3 LO and N 4 LO corrections in the series for 
BjnSR. 



/ 


df 


(df^ECH 


(df l )pMS 


d% st 


(dl st )ECH — c^di 


1 


34.01 


27.25 


28.41 


290 


2561 


2 


26.93 


23.11 


24.09 


203 


1580 


3 


20.21 


19.22 


20.01 


130 


852 


4 


13.84 


15.57 


16.16 


68 


343 


5 


7.83 


12.19 


12.59 


18 


25 


6 


2.17 


9.08 


9.29 


-22 


-130 



Table 4: The results of estimates of the NNLO, N 3 LO and N A LO corrections in the series for 
BjpSR and GLSSR. 
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